An approach to identify gene-environment interactions and reveal new biological insight in complex traits.

There is a long-standing debate about the magnitude of the contribution of gene-environment interactions to phenotypic variations of complex traits owing to the low statistical power and few reported interactions to date. To address this issue, the Gene-Lifestyle Interactions Working Group within the Cohorts for Heart and Aging Research in Genetic Epidemiology Consortium has been spearheading efforts to investigate G × E in large and diverse samples through meta-analysis. Here, we present a powerful new approach to screen for interactions across the genome, an approach that shares substantial similarity to the Mendelian randomization framework. We identify and confirm 5 loci (6 independent signals) interacted with either cigarette smoking or alcohol consumption for serum lipids, and empirically demonstrate that interaction and mediation are the major contributors to genetic effect size heterogeneity across populations. The estimated lower bound of the interaction and environmentally mediated heritability is significant (P < 0.02) for low-density lipoprotein cholesterol and triglycerides in Cross-Population data. Our study improves the understanding of the genetic architecture and environmental contributions to complex traits.

Perceived Stress From Childhood to Adulthood and Cardiometabolic End Points in Young Adulthood: An 18-Year Prospective Study.

BACKGROUND: We investigated how childhood-to-adulthood perceived stress patterns predict adult cardiometabolic risk. METHODS AND RESULTS: This study included 276 participants from the Southern California Children's Health Study (2003-2014), and a follow-up assessment (2018-2021). Perceived stress (Perceived Stress Scale) was initially reported by participants' parents for themselves during early childhood (mean age, 6.3 years), and later self-reported during adolescence (13.3 years) and young adulthood (23.6 years). Participants were grouped into 4 stress patterns: consistently high, decreasing, increasing, and consistently low. Cardiometabolic risk was assessed in young adulthood by carotid artery intima-media thickness, systolic and diastolic blood pressure, obesity, percent body fat, android/gynoid ratio, and glycated hemoglobin. A cardiometabolic risk score was generated by summing the clinically abnormal markers. Multivariable linear and logistic regression models were used to (1) examine the associations between Perceived Stress Scale at 3 time points and adult cardiometabolic risk, and (2) assess the impact of stress pattern on adult cardiometabolic risk. Findings suggested that in adulthood, higher Perceived Stress Scale score was associated with increased overall cardiometabolic risk (ß=0.12 [95% CI, 0.01-0.22]), carotid artery intima-media thickness (ß=0.01 [95% CI, 0.0003-0.02]), systolic blood pressure (ß=1.27 [95% CI, 0.09-2.45]), and diastolic blood pressure (ß=0.94 [95% CI, 0.13-1.75]). Individuals with a consistently high adolescence-to-adulthood stress pattern had greater overall cardiometabolic risk (ß=0.31 [95% CI, 0.02-0.60]), android/gynoid ratio (ß=0.07 [95% CI, 0.02-0.13]), percent body fat (ß=2.59 [95% CI, 0.01-5.17]), and greater odds of obesity (odds ratio, 5.57 [95% CI, 1.62-19.10]) in adulthood, compared with those with a consistently low Perceived Stress Scale score. CONCLUSIONS: Consistently high perceived stress from adolescence to adulthood may contribute to greater cardiometabolic risk in young adulthood.

Genome-Wide Interaction Analysis With DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

BACKGROUND: The Dietary Approaches to Stop Hypertension (DASH) diet score lowers blood pressure (BP). We examined interactions between genotype and the DASH diet score in relation to systolic BP. METHODS: We analyzed up to 9 420 585 single nucleotide polymorphisms in up to 127 282 individuals of 6 population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (n=35 660) and UK Biobank (n=91 622) and performed European population-specific and cross-population meta-analyses. RESULTS: We identified 3 loci in European-specific analyses and an additional 4 loci in cross-population analyses at Pinteraction<5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency, 0.03) and the DASH diet score (Pinteraction=4e-8; P for heterogeneity, 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (Pinteraction=9.4e-7) and 0.20±0.06 mm Hg (Pinteraction=0.001) in Cohorts for Heart and Aging Research in Genomic Epidemiology and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P=4e-273) and cis-DNA methylation quantitative trait loci variants (P=1e-300). Although the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by single nucleotide polymorphisms potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. CONCLUSIONS: We demonstrated gene-DASH diet score interaction effects on systolic BP in several loci. Studies with larger diverse populations are needed to validate our findings.

Trans-ancestral genome-wide association study of longitudinal pubertal height growth and shared heritability with adult health outcomes.

BACKGROUND: Pubertal growth patterns correlate with future health outcomes. However, the genetic mechanisms mediating growth trajectories remain largely unknown. Here, we modeled longitudinal height growth with Super-Imposition by Translation And Rotation (SITAR) growth curve analysis on ~ 56,000 trans-ancestry samples with repeated height measurements from age 5 years to adulthood. We performed genetic analysis on six phenotypes representing the magnitude, timing, and intensity of the pubertal growth spurt. To investigate the lifelong impact of genetic variants associated with pubertal growth trajectories, we performed genetic correlation analyses and phenome-wide association studies in the Penn Medicine BioBank and the UK Biobank. RESULTS: Large-scale growth modeling enables an unprecedented view of adolescent growth across contemporary and 20th-century pediatric cohorts. We identify 26 genome-wide significant loci and leverage trans-ancestry data to perform fine-mapping. Our data reveals genetic relationships between pediatric height growth and health across the life course, with different growth trajectories correlated with different outcomes. For instance, a faster tempo of pubertal growth correlates with higher bone mineral density, HOMA-IR, fasting insulin, type 2 diabetes, and lung cancer, whereas being taller at early puberty, taller across puberty, and having quicker pubertal growth were associated with higher risk for atrial fibrillation. CONCLUSION: We report novel genetic associations with the tempo of pubertal growth and find that genetic determinants of growth are correlated with reproductive, glycemic, respiratory, and cardiac traits in adulthood. These results aid in identifying specific growth trajectories impacting lifelong health and show that there may not be a single "optimal" pubertal growth pattern.

Air pollution and age-dependent changes in emotional behavior across early adolescence in the U.S.

Recent studies have linked air pollution to increased risk for behavioral problems during development, albeit with inconsistent findings. Additional longitudinal studies are needed that consider how emotional behaviors may be affected when exposure coincides with the transition to adolescence - a vulnerable time for developing mental health difficulties. This study investigates if annual average PM2.5 and NO2 exposure at ages 9-10 years moderates age-related changes in internalizing and externalizing behaviors over a 2-year follow-up period in a large, nationwide U.S. sample of participants from the Adolescent Brain Cognitive Development (ABCD) Study®. Air pollution exposure was estimated based on the residential address of each participant using an ensemble-based modeling approach. Caregivers answered questions from the Child Behavior Checklist (CBCL) at the baseline, 1-year follow-up, and 2-year follow-up visits, for a total of 3 waves of data; from the CBCL we obtained scores on internalizing and externalizing problems plus 5 syndrome scales (anxious/depressed, withdrawn/depressed, rule-breaking behavior, aggressive behavior, and attention problems). Zero-inflated negative binomial models were used to examine both the main effect of age as well as the interaction of age with each pollutant on behavior while adjusting for various socioeconomic and demographic characteristics. Against our hypothesis, there was no evidence that greater air pollution exposure was related to more behavioral problems with age over time.

Genome-Wide Interaction Analysis with DASH Diet Score Identified Novel Loci for Systolic Blood Pressure.

Objective: We examined interactions between genotype and a Dietary Approaches to Stop Hypertension (DASH) diet score in relation to systolic blood pressure (SBP). Methods: We analyzed up to 9,420,585 biallelic imputed single nucleotide polymorphisms (SNPs) in up to 127,282 individuals of six population groups (91% of European population) from the Cohorts for Heart and Aging Research in Genomic Epidemiology consortium (CHARGE; n=35,660) and UK Biobank (n=91,622) and performed European population-specific and cross-population meta-analyses. Results: We identified three loci in European-specific analyses and an additional four loci in cross-population analyses at P for interaction < 5e-8. We observed a consistent interaction between rs117878928 at 15q25.1 (minor allele frequency = 0.03) and the DASH diet score (P for interaction = 4e-8; P for heterogeneity = 0.35) in European population, where the interaction effect size was 0.42±0.09 mm Hg (P for interaction = 9.4e-7) and 0.20±0.06 mm Hg (P for interaction = 0.001) in CHARGE and the UK Biobank, respectively. The 1 Mb region surrounding rs117878928 was enriched with cis-expression quantitative trait loci (eQTL) variants (P = 4e-273) and cis-DNA methylation quantitative trait loci (mQTL) variants (P = 1e-300). While the closest gene for rs117878928 is MTHFS, the highest narrow sense heritability accounted by SNPs potentially interacting with the DASH diet score in this locus was for gene ST20 at 15q25.1. Conclusion: We demonstrated gene-DASH diet score interaction effects on SBP in several loci. Studies with larger diverse populations are needed to validate our findings.

Gene-Environment Analyses Reveal Novel Genetic Candidates with Prenatal Tobacco Exposure in Relation to Risk for Childhood Acute Lymphoblastic Leukemia.

BACKGROUND: Associations between maternal tobacco exposure during pregnancy and childhood acute lymphoblastic leukemia (ALL) have yielded mixed results. This may be due to biases in self-reported smoking or other differences in individual-level risk factors. We utilized a biological marker of maternal tobacco exposure to evaluate the association between maternal tobacco exposure during pregnancy, genetics, and subsequent childhood ALL risk in two large population-based studies of childhood ALL in California. METHODS: Maternal exposure to tobacco smoke was assessed with a validated methylation marker (cg05575921) of the aryl hydrocarbon receptor repressor (AHRR) gene in newborn dried blood spots. We adjusted for sex, birthweight, gestational age, mode of delivery, year of birth, AHRR quantitative trait locus (mQTL) rs77111113, and a polygenetic risk score for childhood ALL. We additionally adjusted for principal components in a gene-environment interaction testing method that incorporates gene-only and environment-only effects along with interactions. RESULTS: AHRR hypomethylation overall was not associated with childhood ALL. In gene-environment interaction testing, several genetic variants displayed significant interaction with AHRR hypomethylation and childhood ALL. CONCLUSIONS: Our results suggest that novel candidates in PTPRK and DPP6 may play a role in tobacco-related leukemogenesis. Further research is necessary to better understand the effects of tobacco and these variants on childhood ALL risk. IMPACT: Despite the lack of an overall "main effect," tobacco exposure during pregnancy affects childhood ALL risk depending on specific genetic variants.

Sex-specific effects in how childhood exposures to multiple ambient air pollutants affect white matter microstructure development across early adolescence.

Ambient air pollution is ubiquitous, yet questions remain as to how it might impact the developing brain. Large changes occur in the brain's white matter (WM) microstructure across adolescence, with noticeable differences in WM integrity in male and female youth. Here we report sex-stratified effects of fine particulate matter (PM2.5), nitrogen dioxide (NO2), and ozone (O3) on longitudinal patterns of WM microstructure from 9-13 years-old in 8,182 (49% female) participants using restriction spectrum imaging. After adjusting for key sociodemographic factors, multi-pollutant, sex-stratified models showed that one-year annual exposure to PM2.5 and NO2 was associated with higher, while O3 was associated with lower, intracellular diffusion at age 9. All three pollutants also affected trajectories of WM maturation from 9-13 years-old, with some sex-specific differences in the number and anatomical locations of tracts showing altered trajectories of intracellular diffusion. Concentrations were well-below current U.S. standards, suggesting exposure to these criteria pollutants during adolescence may have long-term consequences on brain development.

A note on p-value multiple-testing adjustment for two-step genome-wide gene-environment interactions scans.

Two-step testing is the state-of-the art approach for performing genome-wide interaction scans (GWIS). It is computationally efficient and yields higher power than standard single-step-based GWIS for virtually all biologically plausible scenarios. However, while two-step tests control the genome-wide type I error rate at the desired level, the lack of associated valid p-values can make it difficult for users to compare with single step-results. We show how multiple-testing adjusted p-values can be defined for two-step test based on standard multiple-testing theory, and how they can be in turn scaled to make valid comparisons with single-step tests possible.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Air pollution and emotional behavior in adolescents across the U.S.

Recent studies have linked air pollution to increased risk for behavioral problems during development, albeit with inconsistent findings. Additional longitudinal studies are needed that consider how emotional behaviors may be affected when exposure coincides with the transition to adolescence - a vulnerable time for developing mental health difficulties. This study examines how annual average PM2.5 and NO2 exposure at ages 9-10 years relates to internalizing and externalizing behaviors over a 2-year follow-up period in a large, nationwide U.S. sample of participants from the Adolescent Brain Cognitive Development (ABCD) Study®. Air pollution exposure was estimated based on the residential address of each participant using an ensemble-based modeling approach. Caregivers answered questions from the Child Behavior Checklist (CBCL) at baseline and annually for two follow-up sessions for a total of 3 waves of data; from the CBCL we obtained scores on internalizing and externalizing problems plus 5 syndrome scales (anxious/depressed, withdrawn/depressed, rule-breaking behavior, aggressive behavior, and attention problems). Zero-inflated negative binomial models were used to examine both the main effect of age as well as the interaction of age with each pollutant on behavior while adjusting for various socioeconomic and demographic characteristics. Overall, the pollution effects moderated the main effects of age with higher levels of PM2.5 and NO2 leading to an even greater likelihood of having no behavioral problems (i.e., score of zero) with age over time, as well as fewer problems when problems are present as the child ages. Albeit this was on the order equal to or less than a 1-point change. Thus, one year of annual exposure at 9-10 years is linked with very small change in emotional behaviors in early adolescence, which may be of little clinical relevance.

Longitudinal hierarchical Bayesian models of covariate effects on airway and alveolar nitric oxide.

Biomarkers such as exhaled nitric oxide (FeNO), a marker of airway inflammation, have applications in the study of chronic respiratory disease where longitudinal studies of within-participant changes in the biomarker are particularly relevant. A cutting-edge approach to assessing FeNO, called multiple flow FeNO, repeatedly assesses FeNO across a range of expiratory flow rates at a single visit and combines these data with a deterministic model of lower respiratory tract NO to estimate parameters quantifying airway wall and alveolar NO sources. Previous methodological work for multiple flow FeNO has focused on methods for data from a single participant or from cross-sectional studies. Performance of existing ad hoc two-stage methods for longitudinal multiple flow FeNO in cohort or panel studies has not been evaluated. In this paper, we present a novel longitudinal extension to a unified hierarchical Bayesian (L_U_HB) model relating longitudinally assessed multiple flow FeNO to covariates. In several simulation study scenarios, we compare the L_U_HB method to other unified and two-stage frequentist methods. In general, L_U_HB produced unbiased estimates, had good power, and its performance was not sensitive to the magnitude of the association with a covariate and correlations between NO parameters. In an application relating height to longitudinal multiple flow FeNO in schoolchildren without asthma, unified analysis methods estimated positive, statistically significant associations of height with airway and alveolar NO concentrations and negative associations with airway wall diffusivity while estimates from two-stage methods were smaller in magnitude and sometimes non-significant.

Ambient temperature and air pollution associations with suicide and homicide mortality in California: A statewide case-crossover study.

BACKGROUND: Higher ambient temperature and air pollution may contribute to increased risk of behaviors harmful to oneself or to others; however, quantitative evidence is limited. We examined the relationship of deaths due to suicide and homicide with temperature and air pollution in California-a state prone to high levels of both exposures. METHOD: California death certificates from 2014 to 2019 were used to identify deaths due to suicide and homicide. Residential data for decedents were used to assign exposure to daily temperature (maximum[Tmax], minimum[Tmin]) and daily average air pollution concentrations (particulate matter <10 µm[PM10] and < 2.5 µm[PM2.5], nitrogen dioxide[NO2], ozone[O3]). Tmin served as a surrogate for nighttime temperature. A time-stratified case-crossover study design using conditional logistic regression was used to assess the effects of daily exposure to temperature and air pollutants on suicide and homicide mortality, adjusting for relative humidity. Effect modification by sex and age was assessed. RESULTS: We observed 24,387 deaths due to suicide and 10,767 deaths due to homicide. We found a monotonic temperature association for both outcomes. A 5 °C increase in Tmax at lag-2 and Tmin at lag-0 was associated with 3.1 % (95 % confidence interval [CI]: 1.1 %-5.2 %) and 3.8 % (95%CI: 0.9 %-6.8 %) increased odds of death due to suicide, respectively. The increased odds of homicide mortality per 5 °C increase in Tmax at lag-0 and Tmin at lag-1 were 4.9 % (95%CI: 1.6 %-8.1 %) and 6.2 % (95%CI: 1.6 %-11.0 %), respectively. No air pollutant associations were statistically significant. Temperature associations were robust after adjustment for PM2.5. Some temperature effects were larger among women for suicide and men for homicide mortality, and among those over age 65 years for both outcomes. CONCLUSION: Risk of suicide and homicide mortality increases with increasing daily ambient temperatures. Findings have public health relevance given anticipated increases in temperatures due to global climate change.

Improved two-step testing of genome-wide gene-environment interactions.

Two-step tests for gene-environment ( G × E $G\times E$ ) interactions exploit marginal single-nucleotide polymorphism (SNP) effects to improve the power of a genome-wide interaction scan. They combine a screening step based on marginal effects used to "bin" SNPs for weighted hypothesis testing in the second step to deliver greater power over single-step tests while preserving the genome-wide Type I error. However, the presence of many SNPs with detectable marginal effects on the trait of interest can reduce power by "displacing" true interactions with weaker marginal effects and by adding to the number of tests that need to be corrected for multiple testing. We introduce a new significance-based allocation into bins for Step-2 G × E $G\times E$ testing that overcomes the displacement issue and propose a computationally efficient approach to account for multiple testing within bins. Simulation results demonstrate that these simple improvements can provide substantially greater power than current methods under several scenarios. An application to a multistudy collaboration for understanding colorectal cancer reveals a G × Sex interaction located near the SMAD7 gene.

Differences in Self-Reported Food Allergy and Food-Associated Anaphylaxis by Race and Ethnicity Among SAPPHIRE Cohort Participants.

BACKGROUND: Although food allergies are considered common, relatively little is known about disparities in food allergy by race in the United States. OBJECTIVE: To evaluate differences in reported food allergy and food-associated anaphylaxis among individuals enrolled in a longitudinal cohort study from metropolitan Detroit, Michigan. METHODS: Participants in the Study of Asthma Phenotypes and Pharmacogenomic Interactions by Race-Ethnicity (SAPPHIRE) were asked about food allergies, including the inciting food and associated symptoms. Individuals were considered to have food-associated anaphylaxis if symptoms coincided with established clinical criteria. Logistic regression was used to assess whether race difference persisted after adjusting for and stratifying by potential confounders. African genetic ancestry was individually estimated among African American SAPPHIRE participants to assess whether ancestry was associated with food allergy. RESULTS: Within the SAPPHIRE cohort, African American participants were significantly more likely to report food allergy (26.1% vs 17%; P = 3.47 × 10-18) and have food-associated anaphylactic symptoms (12.7% vs 7%; P = 4.65 × 10-14) when compared with European American participants. Allergy to seafood accounted for the largest difference (13.1% vs 4.6%; P = 1.38 × 10-31). Differences in food allergy by race persisted after adjusting for potential confounders including asthma status. Among African American participants, the proportion of African ancestry was not associated with any outcome evaluated. CONCLUSION: Compared with European Americans, African Americans appear to be at higher risk for developing food allergy and food-associated anaphylaxis, particularly with regard to seafood allergy. The lack of association with genetic ancestry suggests that socioenvironmental determinants may play a role in these disparities.

Multi-omics insights into the biological mechanisms underlying statistical gene-by-lifestyle interactions with smoking and alcohol consumption.

Though both genetic and lifestyle factors are known to influence cardiometabolic outcomes, less attention has been given to whether lifestyle exposures can alter the association between a genetic variant and these outcomes. The Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium's Gene-Lifestyle Interactions Working Group has recently published investigations of genome-wide gene-environment interactions in large multi-ancestry meta-analyses with a focus on cigarette smoking and alcohol consumption as lifestyle factors and blood pressure and serum lipids as outcomes. Further description of the biological mechanisms underlying these statistical interactions would represent a significant advance in our understanding of gene-environment interactions, yet accessing and harmonizing individual-level genetic and 'omics data is challenging. Here, we demonstrate the coordinated use of summary-level data for gene-lifestyle interaction associations on up to 600,000 individuals, differential methylation data, and gene expression data for the characterization and prioritization of loci for future follow-up analyses. Using this approach, we identify 48 genes for which there are multiple sources of functional support for the identified gene-lifestyle interaction. We also identified five genes for which differential expression was observed by the same lifestyle factor for which a gene-lifestyle interaction was found. For instance, in gene-lifestyle interaction analysis, the T allele of rs6490056 (ALDH2) was associated with higher systolic blood pressure, and a larger effect was observed in smokers compared to non-smokers. In gene expression studies, this allele is associated with decreased expression of ALDH2, which is part of a major oxidative pathway. Other results show increased expression of ALDH2 among smokers. Oxidative stress is known to contribute to worsening blood pressure. Together these data support the hypothesis that rs6490056 reduces expression of ALDH2, which raises oxidative stress, leading to an increase in blood pressure, with a stronger effect among smokers, in whom the burden of oxidative stress is greater. Other genes for which the aggregation of data types suggest a potential mechanism include: GCNT4×current smoking (HDL), PTPRZ1×ever-smoking (HDL), SYN2×current smoking (pulse pressure), and TMEM116×ever-smoking (mean arterial pressure). This work demonstrates the utility of careful curation of summary-level data from a variety of sources to prioritize gene-lifestyle interaction loci for follow-up analyses.

Interactions between folate intake and genetic predictors of gene expression levels associated with colorectal cancer risk.

Observational studies have shown higher folate consumption to be associated with lower risk of colorectal cancer (CRC). Understanding whether and how genetic risk factors interact with folate could further elucidate the underlying mechanism. Aggregating functionally relevant genetic variants in set-based variant testing has higher power to detect gene-environment (G × E) interactions and may provide information on the underlying biological pathway. We investigated interactions between folate consumption and predicted gene expression on colorectal cancer risk across the genome. We used variant weights from the PrediXcan models of colon tissue-specific gene expression as a priori variant information for a set-based G × E approach. We harmonized total folate intake (mcg/day) based on dietary intake and supplemental use across cohort and case-control studies and calculated sex and study specific quantiles. Analyses were performed using a mixed effects score tests for interactions between folate and genetically predicted expression of 4839 genes with available genetically predicted expression. We pooled results across 23 studies for a total of 13,498 cases with colorectal tumors and 13,918 controls of European ancestry. We used a false discovery rate of 0.2 to identify genes with suggestive evidence of an interaction. We found suggestive evidence of interaction with folate intake on CRC risk for genes including glutathione S-Transferase Alpha 1 (GSTA1; p = 4.3E-4), Tonsuko Like, DNA Repair Protein (TONSL; p = 4.3E-4), and Aspartylglucosaminidase (AGA: p = 4.5E-4). We identified three genes involved in preventing or repairing DNA damage that may interact with folate consumption to alter CRC risk. Glutathione is an antioxidant, preventing cellular damage and is a downstream metabolite of homocysteine and metabolized by GSTA1. TONSL is part of a complex that functions in the recovery of double strand breaks and AGA plays a role in lysosomal breakdown of glycoprotein.

Transcriptomic Response to Calcium in Normal Colon Organoids is Impacted by Colon Location and Sex.

Observational studies indicate that calcium supplementation may protect against colorectal cancer. Stratified analyses suggest that this protective effect may differ based on anatomic subsite and sex, but these hypotheses have been difficult to test experimentally. Here, we exposed 36 patient-derived organoid lines derived from normal colon biopsies (21 right colons, 15 left colons) of unrelated subjects (18 female, 18 male) to moderate (1.66 mmol/L) or high (5.0 mmol/L) concentrations of calcium for 72 hours. We performed bulk RNA-sequencing to measure gene expression, and cell composition was inferred using single-cell deconvolution in CIBERSORTx. We tested for significant differences in gene expression using generalized linear models in DESeq2. Exposure to higher levels of calcium was associated with changes in cell composition (P < 0.05), most notably increased goblet and reduced stem cell populations, and differential expression of 485 genes (FDR < 0.05). We found that 40 of these differentially expressed genes mapped to genomic loci identified through colorectal cancer genome-wide association studies, suggesting a potential biologic overlap between calcium supplementation and inherited colorectal cancer risk. Stratified analyses identified more differentially expressed genes in colon organoids derived from right sided colon and male subjects than those derived from left sided colon and female subjects. We confirmed the presence of a stronger right-sided effect for one of these genes, HSD17B2 using qPCR in a subset of matched right and left colon organoids (n = 4). By relating our findings to genetic data, we provide new insights into how nutritional and genetic factors may interact to influence colorectal cancer risk. PREVENTION RELEVANCE: A chemopreventive role for calcium in colorectal cancer is still unclear. Here, we identify mechanisms through which calcium supplementation may reduce risk. Calcium supplementation increased differentiation and altered expression of colorectal cancer-related genes in a large study of patient-derived colon organoids. These findings were influenced by colon location and sex.

Gene-lifestyle interactions in the genomics of human complex traits.

The role and biological significance of gene-environment interactions in human traits and diseases remain poorly understood. To address these questions, the CHARGE Gene-Lifestyle Interactions Working Group conducted series of genome-wide interaction studies (GWIS) involving up to 610,475 individuals across four ancestries for three lipids and four blood pressure traits, while accounting for interaction effects with drinking and smoking exposures. Here we used GWIS summary statistics from these studies to decipher potential differences in genetic associations and G×E interactions across phenotype-exposure-ancestry combinations, and to derive insights on the potential mechanistic underlying G×E through in-silico functional analyses. Our analyses show first that interaction effects likely contribute to the commonly reported ancestry-specific genetic effect in complex traits, and second, that some phenotype-exposures pairs are more likely to benefit from a greater detection power when accounting for interactions. It also highlighted modest correlation between marginal and interaction effects, providing material for future methodological development and biological discussions. We also estimated contributions to phenotypic variance, including in particular the genetic heritability conditional on the exposure, and heritability partitioned across a range of functional annotations and cell types. In these analyses, we found multiple instances of potential heterogeneity of functional partitions between exposed and unexposed individuals, providing new evidence for likely exposure-specific genetic pathways. Finally, along this work, we identified potential biases in methods used to jointly meta-analyze genetic and interaction effects. We performed simulations to characterize these limitations and to provide the community with guidelines for future G×E studies.